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A NEW PROOF OF SEMICIRCLE LAW OF FIXED TRACE 
SQUARE ENSEMBLE 

DA XU AND LIHE WANG 



Abstract. In the present paper, we give a simple proof of the level density of 
fixed trace square ensemble. We derive the integral equation of the level density 
of fixed trace square ensemble. Then we analyze the asymptotic behavior of the 
5H level density. 

in 

,__! 1. Introduction 

,_i E.P.Wigner in [TT] and [H] showed that the level density (one point correlation) 

rS^ of "border matrix" ensemble which is a kind of real symmetric matrix ensemble, 

Mh asymptotically approaches to the semicircle, i.e., 

-S (1-1) '^sys,n{x) -^ ^/2N-xyTT. 

^ Porter and Rosenzweig also found that the level density of some special random 

G matrices satisfy semi-circle law [7j |8] . 

Wigner's movitation arose from the consideration of properties of wave functions 
^SJ of quantum mechanical system that are extremely complicated that statistical con- 

^ siderations can be applied to them. From a computational point of view, random 

matrix theory, as a statistical theory, is very important to describe extremely com- 
^^ plicated quantum mechanical systems such as large atoms, QCD physics(see [TD] 

fvj for review) . It is instructive to continue to work on the properties of RMT in the 

problems of complicated quantum systems. 

Recently there has been work on the analysis of top eigenvalues of random 
QQ matrices ([13]). When we evaluate some quantities, for example the ratio of the ex- 

f^ pectation of the largest square of eigenvalues and the expectation of the square of 

an arbitrary eigenvalues 

^ {ma.xi<i<Nxf)GUE _ /• ••/maxi<i<Ara;fe ^i ^^ Yl^^j{xi - Xj)^dxidx2- ■ ■ dxN 

C^ {xi)gue ^ ■ ■■ ^ Tiia.yii<i<Nx\e'^i ■^i]\^^.{x^- XjYdxidx2- --dxM 

, we can change the integrals into polar coordinates. We can see the integration 
on vandermonde space is the difficulty. Therefore it is important to consider the 
"fixed trace square ensemble" on M^ 

on the A^— 1 dimensional unit sphere {(a;i, • • • , a; at); ^j^ xf — 1} and where C^ is a 
normalization constant. Stieltjes in 1914 proved a theorem that states the maximum 
value of vandermonde polynomial J^,- \xi — Xj\'^ in the ball {x\ + x"^- ■ ■ -\~ x'f^ < 



(1.3) Pv,n{xi,X2,-- ■ ,xn) = CvY\_\ 
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N{N—1)} achieves its maximal values at zeros of Hermite polynomials. In [13], the 
author considered the bounded trace square ensemble and obtained semicircle law 
of bounded trace square ensemble. In [?], the authors reobtained the semicircle law 
of fixed trace square ensemble by comparing it with GUE and saddle point method. 
In the present paper, we would like to give a simple proof of the semicircle law of 
fixed square ensemble. 

2. INTEGRALS OF FIXED TRACE SQUARE ENSEMBLE 

Let's give a general definition of the ensemble on a homogeneous space. 

Definition 1. P{X) — P{xi, X2,- ■ ■ , x^) is a probability distribution satisfying on 
a homogenous space M with respect to the permutation group action on xi, X2, • • • , xn , 
which are the local coordinates. The level density is defined to be 



(2.1) ap{x) = N Jp{xi,X2,- 



,Xi\[)dx2 ■ ■ ■ dx 



N- 



In this paper, we are concerned with the case that M is the A^ — 1 dimensional 
sphere. 

We would like to compute the level density of the ensemble P{xi, X2, • • • , xn) = 
Cjv ni<j \^i~^j\^ oil the N dimensional ball Bi — {(xi, ^2, • • • ,xn); J2i=i ^f < !}• 
The level density on sphere will easily follow. We will use Selberg's integral [9]. 

Theorem 2.1. (Selberg) For any positive integer n, let dx = dxi ■ ■ ■ dxN, 
A{x) = A{xi,---,X2)= Y\_ (^^i -2;,), i/ n > 1, 

l<j<l<n 

(2.2) A(x) = l, ifn^l, 
and 

N 

(2.3) <Pix) EE <P{x^, ■■■,xn) = |A(x)|2^ n ^r'(l - ^^)^''- 

i=i 
Then 

(2.4) /(a, I3,j,n) = J^ ■■■ j^ ^{x)dx =\{f^ 



" ' r( l+7 + J7)r(aj-7)r(/? + J7 
^.^0 - .- + 7)r(a + P+{n + j- 1)7 ' 



and for I < m < n, 

(2.5) 

xiX2 ■ ■ ■ x^<P{x)dx = TT " ;" — ^'."^ / •••/ <P{x)dx, 

valid for integer n and complex a,(3,"f with 

1 Bea BeB 

(2.6) i?ea > 0, i?e/3 > 0, i?e7 > -min(-, -, ^). 

n n — 1 n— 1 

We need to compute the constant in the expression of the ensemble. We cannot 
directly apply Selberg's theorem, because the domain we integrate on is the unit 
ball instead of M.^ . By the virtue of gamma function, we can change the integral 
on the unit ball to the integral on M^ of Selberg's type. Let us prove the following 
identity first. 
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Proposition 1. For any R> 0, we have 



.7) / A'dx = R-' ^} (2.)"/^2-^V2 J^ r(l + j) 



Proof. In Selberg's integral, letting Xi — ^ + ^p] and L — > oo,then 



2L '^ 2l 

N N 



(2.8) r ■■■ r |A|2^ n e-'^^Ux^ = (2^)"/2(2a)-^(^("-i)+i)/2 f[ ^^^ + ^'^^ 



Multiplying a*^ ^e ''to (2.8 land integrating with respect to a, we get 

N N 



-OO >/— OO ■_i ■_i 



JV 



r(/3 - iV(7(n - 1) + l)/2) . .„/2.-^(.(n-i)+i)/2 TT r(l+J7) 

^ ■ ^ r(/3) ^ ^ IJ r(i + 7)- 

We change variables by x, — , ^' ^^. The determinant is 

n-lj 1 1 1 

(2.11) -^ cir, _ i i _ i 



ry-'dry (1 _ ^f^^ y^^)-^ (1 - ^f^^ j;|)f (1 - ^f^^ j;|)t+i ' 

where rx,ry are the radius coordinates in the polar coordinate system. Therefore 
we have 



(2.12) / .../ \A\'-^{l + Y.<r'X{dx, 

J —OO J —OO ■_■, ■_^ 



i=i j=i 



OO J — OO _ 

A^ A^ 



ii) 



(2.13) = / |A(y)p'^(i-^y,2)^-^W"~i)+i)/2)-in^y 

Then we end up with 

/ |A(y)p'^(l - J2 y^r-^M-^^+^y^)-^ fj dy, 

_ r(/3 - iV(7(n - 1) + l)/2) ,, ^^(,(„_,,^,,,, ^ r(^ 

Now set 7 = 1 and /? - N{-f{n - 1) + l)/2) -1 = 0. We get 

N N 

By rescaling, we conclude the result of this proposition. D 

we are now ready to prove the following identity. 
Theorem 2.2. The level density of P{x) = const. A^(x) in the unit hall is 
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Proof. The level density of P{x) is 

CTp(2/i) =N f ■■■ f dy2--- dvNA^iy) 

J Jvi+vi+---+y%<'i--y^ 

= N ■■■ dy2--- dyN ^(yi - y^f J| (?/» - Vj) 



yi+yi+---+y%<i-y^ j=2 2<i<j<N 

= Nf---f dy2---dyNf[{yl-2yiyj+y]) [] {y^ - yjf 

•I •'yl+yl+-+yl<'^'y'' j=2 2<i<j<N 



'yi+yi + ---+yN<l-y^ p=0 g=0 2<i<i<iV 

(2.16) 

Af-lTV-l-p 
p=0 ?=0 

D 



One can easily apply Aomoto's technique ([T]) to compute (xiX2 • • • Xm)- How- 
ever, we wish that we could get an analytic expression of the level density. Unfortu- 
nately, as wee have seen the above theorem, we have to explicitly write expectations 
such as {x\x'^ ■ ■ ■ x'^^Xm+iXm+2 ' ' ' Xn) which turns out to be very hard and remains 
an open problem. This fact makes us have to think about the numerical solution 
of the level density. 



3. THE INTEGRAL EQUATION OF LEVEL DENSITY 

We would like to derive the integral equation of the level density of fixed trace 
square ensemble in another way. The idea is this. Consider the level density 
of GUE. On the sphere with radius r, the level density is -.^{-) with weight 
Cj^^e-''V"-V"("-i) where 

(3.1) Cj, = r e-V^- V^(^-i)dr = ir(^). 



So we have 

Theorem 3.1. The level density of fixed trace square ensemble satisfies the integral 
equation 

(3.2) 2r-i( — )/ e-^-r^ -^a^,^{x/r)dr = aGUE.n{^)- 
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Proof. For any f{xi) £ C^{R), we change the following integral into polar coordi- 
nate sytem. 



(^GUEixi)f{xi)dxi 
N / CGUE,Ne^ ^1 ''' Wi^i - XjYf{xi)dxidx2 ■ ■ -dx 



--N / e-''V^-i+^(^-i)dr / ^\n)f{rn^)dn 

Jo Jn 

^OO pi 





oo 



const. / e-'' r^-i+^(^-i)-idr / dx^a,{xjr)f{x{) 







(3.3) =2r-i / e-'-V^-i+w(^-i)-idr / dxia^{xi/r)f{xi). 



The constant in (3.3) is computed by Proposition 1. D 

Similarly, we can derive the integral equation of two point cluster function of 
fixed trace square ensemble Yy^2{xi,X2) in terms of the cluster function of GUE 

1^ / \ / sinTrfxi— 3:2) \2 

yGUE,2{^l,X2) = ( ^(i,-^,)0 ■ 

It immediately follows that the asymptotic behavior of (7^,(0). 
Corollary 1. 

(3.4) -^-(«)-^f^^(l + «W)- 

4. MAIN THEOREM 

It is well known that the eigenvalue distribution of GUE, p{x) — ^/ j^aGUE,NW'2Nx) 



N 

goes to semicircle: 

/A ^\ / \ n/ \ f \/l — x'^, if \x\ < 1 

(4.1) p{x) -. Six) = I V ^^ ^yj^ - 

Theorem 4.1. It is possible to prove the asymptotic behavior of the integral equa- 
tion (3.2). In fact, The normalized eigenvalue distribu tion of fixed trace square 
ensemble Oy m{x) = -^(Jv n{-^) Qoes to the semicircle (A.l). 

Proof. First we have to estimate the L°° norm of ctjv. 
Define a function 

/■oo 

(4.2) a^)^mma,.M{-)2V-'{N^l2) e'^\''' -'' dr} 

t. J\x\ 

Note that £,{x) is a continuous function on [0, 00). Then 77^:^ (at x — 0, this function 
is defined to be zero) is also a continuous function. Since (7gue,n — > 0(|1|) as 
x —^ 00, then 

(A Q^ ( let \\ CfGUE,N{x) 

(4.3) cr„,Ar(x/^x)) - • " 



2r-i(7V2/2) /|^ e-'-V^'-2dr 
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as X ^ oo Therefore cr^j^^-^rp-) — > and j^ —>■ 1. By the expansion of herniite 
functions in |3], we know Vx £ [0, 1], 



(4 4) ^--^^^^ < CN'2 



Therefore 

(4.5) cTy^N<CN^. 

We shall use the following lemma: 
Lemma 4.2. VO < a < 1, V?n > 0, if N is sufficiently large, 



(4.6) 1 - N- < -Voo--;.^.., dr < 1. 



J M^ 

rOO 

Jo e- 



lN/V2-Na "^ ^ 



The proof is straightforward. 

For \/a > 0, if A^ is sufficiently large, 



2r-^(iVV2) / e-"- r" -'a, N{x/r)dr 

J\x\ 

=2r-i(iVV2)(/ +/ +/ )e-'-V^'-V,,^(x/r)dr 

J\x\ JN/^/2'-Na JN/^/2+Na 

pN/V^-Na poo pN/V^+Na 

=2r-i(iVV2)(/ +/ )e--\''"-^N'-^dr+ e-'-'r^'" Vw(a;/r)dr 

J\x\ JN/\/2+Na JN/^/2-Na 



^2T-'{N'/2) / e-"- r^^ -V,,jv(^A)dr + 0(1/7V") 

(4.7) 

-2r-i(7VV2)a„,w(x/r(x)) / e-'^V^'-2^r + 0(l/iV5), 



where iV/v/2 - Na< r{x) < N/V^ + Na. Therefore VO < y, 



(4.8) a, Ay) = ^, Zl!^^''^^ <yGUE,N{r{y)y) + o{i/n^). 



We notice that ,^^^^2^S ,. ~^ N/V2 + 0(l/y/N)(see [5"). Therefore, 
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\Pv.N{y)- S{y)\ 

TV 2 VjV 



<\^^<JGUE,N{r2y/VN) - Siy)\ + o(l) 

<\\l ^<^GUEM^^ir^y/i^^^))) - S{P2y/{VNV2N))\ + \S{r2y/{VNV2N)) - S{y)\ + o(l) 

(4.9) 

< sup \S{x)^S{y)\ + o{l) 

\x~y\<a 

D 

We use the software Matlab to estimate the level density ct^^at of fixed trace 
square ensemble in the case A^ = 10, A^ = 50, and N = 100. Surprisingly, these 
graphs look really like semicircles. 



DA XU AND LIKE WANG 




-1 -0.8 -O.G -0.4 -0.2 0.2 0.4 O.G 0.8 1 



Figure 1. Level density of fixed trace square ensemble N=10 




Figure 2. Level density of fixed trace square ensemble N=:50 
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Figure 3. Level density of fixed trace square ensemble N=100 
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Figure 4. the coincidence of the semicircle and the level density 
when iV == 100 
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